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© METHODE POUR CONSTITUER UN MODELE REPRESENTATIF D'ECOULEMENTS POLYPHASIQUES DANS 
DES CONDUITES DE PRODUCTION PETROLIERE. 

Q>7) La methode proposee a pour objet de constituer un 
modele repr6sentatif des ecoulements permanents et tran- 
sitoires dans des conduites d'un melange de fluides poly- 
phasiques, qui prend en compte un ensemble de variables 
definissant les proprietes des fluides et des modes d'ecou- 
lement: a phases separees, disperses, intermittents, et 
aussi des dimensions et pentes des conduites d'achemine- 
ment. On determine les grandeurs caracterisant recoup- 
ment en resolvant un ensemble d'equations de transport, 
une equation de conservation de masse par constituant et 
une equation de quantite de mouvement du melange, et en 
utilisant un module hydrodynamique et une caracteristique 
thermodynamique des fluides. On forme le modele en 
considerant que le melange est sensiblement a I'equilibre a 
chaque instant et que la composition du melange polypha- 
sique est variable tout le long de la conduite. 

Application a I'etude de reseau d'acheminement d'hydro- 
carbures par exemple. 
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La presente invention concerne une methode pour constituer un modele 
representatif des ecoulements stationnaires et transitoires dans des conduites d'un 
melange de fluides polyphasiques. 

5 La moderation realisee par la methode selon l'invention, permet de tenir 

compte des phenomenes de transfert de masse entre phases et de quantites de 
mouvement entre les phases du melange, a partir d'un ensemble de variables 

definissant— les— proprietES-desHluidesr^eure 

variations de la pente des conduites par rapport a 1'horizontale. Le modele ainsi 

10 constitue permet de faciliter la conception de reseaux de transfert d'effluents 
petroliers par exemple. 

La methode selon l'invention convient notamment pour modeliser le 
comportement de melanges polyphasiques d'hydrocarbures circulant dans des 
pipelines depuis des sites d'exploitation de gisements jusqu'a des sites de chargement 
15 ou de traitement par exemple. 

II est bien connu des specialistes que les modes d'ecoulement de fluides 
polyphasiques dans des tubes sont extremement varies et complexes. Les 
ecoulements diphasiques, par exemple, peuvent etre stratifies, la phase liquide 
s'ecoulant dans la partie inferieure de la conduite, ou intermittents avec une 

20 succession de bouchons liquides et gazeux, ou bien encore disperses, le liquide etant 
entraine sous forme de fines gouttelettes. Le mode d'ecoulement varie notamment 
avec Tinclinaison des conduites par rapport a l'horizontale et il depend du debit de la 
phase gazeuse, de la temperature etc. Le glissement entre les phases qui varie selon 
que Ton considere les tron?ons ascendants ou descendants, entraine des variations de 

25 pression sans qu'il y ait toujours compensation. Les caracteristiques du reseau 
d'ecoulement (dimensionnement, pression, debit de gaz etc.) doivent etre 
determinees avec soin. 

Parmi les nombreuses publications portant sur le comportement des 
ecoulements notamment diphasiques dans des conduites, on peut citer par exemple : 
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- Fabre, J., et al 1983, Intermittent gas-liquid flow in horizontal or slightly inclined 
pipes, Int. Conference on the Physical Modelling of Multi-Phase Flow, Coventry, 
England, pp 233, 254; ou 

- Fabre, J., et al 1989, Two fluid/two flow pattern model for transient gas liquid 
5 flow in pipes, Int. Conference on Multi-Phase Flow, Nice, France, pp 269, 284, 

Cranfield, BHRA. 

Une methode de modelisation existe ou Ton traite les changements de phase 
par des processus iteratifs : Tetat du melange est cense connu a priori et, si~cela~ 
conduit, a des incoherences, apres calcul hydrodynamique, les calculs sont repetes 
10 avec un nouvel etat du melange. Cette methode est lourde et peut etre la source de 
problemes de convergence. 

Une methode de modelisation appliquee a des milieux poreux est decrite par 
exemple par Eymard R., Gallouet T., 1991, Traitement des changements de phases 
dans la modelisation de gisements petroliers. Journees numeriques de Besangon, 23 
15 et 24 septembre 1991. 

Par le brevet US 5 550 761 du demandeur, on connait une methode de 
modelisation des 6couIernents polyphasiques permanerits ou transitoires qui prend en 
compte un ensemble de variables definissant les proprietes des fluides et des modes 
d'ecoulement, et aussi des dimensions et pentes des conduites d'acheminement. On 
20 determine les grandeurs caracterisant l'ecoulement en resolvant un ensemble 
d'equations de transport avec une equation de conservation par phase et une equation 
de quantite de mouvement du melange, et en utilisant un modele hydrodynamique et 
une caracteristique thermodynamique des fluides. 

Pour obtenir ce modele hydrodynamique, on caracterise les regimes 
25 d'ecoulement par un parametre variant entre 0 et 1 representatif de la fraction de 
Tecoulement qui est dans un etat separe (les phases sont stratifiees verticalement ou 
radialement par exemple), on determine tout regime d'ecoulement au cours de la 
resolution des equations de transport, par une comparison de la valeur courante de 
la fraction liquide dans les bouchons et de celle des regions a mode d'ecoulement 
30 disperse, on determine aussi la vitesse des bouchons de la phase gazeuse par rapport 
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a une vitesse critique, et on impose durant la resolution des relations de fermeture, 
des contraintes de continuite aux frontieres entre les regimes, aux fractions de 
volume gazeux ainsi qu ! a la vitesse de deplacement des bouchons. 

Dans la methode precedents on a choisi une approche "par phase" (liquide- 
5 gaz) ou la conservation de la masse est traduite par une equation de conservation par 
phase et le transfert de masse entre phases est exprime par un terme de desequilibre 
proportionnel a l'ecart entre deux valeurs de fraction massique vapeur, Tune jmva e(l , 

correspondant a I'equilibre, qui est fournie par la thermodynamique a composition 
gl nhale constante, l'autre etant celle calculee en tenant compte du g lissement entre 
10 phases : 



ou AGKL est un facteur dependant a priori du fluide et de la configuration de 
recoupment. 

On a constate a Tusage qu'il est difficile de definir une formulation de ce 
15 terme de desequilibre, qui s'applique a toutes les situations: pentes locales des 
conduites avec des points hauts et des points bas, importants transferts de masse 
entre phases. II n'existe pas de methode fiable et robuste permettant de prendre 
correctement en compte le terme de desequilibre entre les phases; 1'approche liquide- 
gaz ou "par phase" ne donne pas de resultats satisfaisants pour traiter les cas ou se 
20 produisent d'importants transferts entre phases. 

Pour constituer un modele representatif des ecoulements permanents et 
transitoires dans des conduites d'un melange polyphasique, qui prend en compte un 
ensemble de variables definissant les proprietes des fluides et des modes 
d'6coulenient: a phases separees, disperses, intermittents, et aussi des dimensions et 
25 pentes des conduites d'acheminement, la methode selon Tinvention comporte 
1'utilisation d f un modele hydrodynamique de type a flux de derive et d'un modele 
thermodynamique pour definir les proprietes des constituants et la resolution d'un 
ensemble d'equations de conservation de la masse par constituant, de conservation de 
la quantite de mouvement du melange et de transfert d'energie dans le melange. 



m G = AKGL- fmva eq - 
V 



Pa'Xr.-Vc ^ 

Pg'RcVg+Pl'Rl'VJ 
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La methode selon l'invention est caracterisee en ce que Ton forme le modele 
en considerant que le melange est sensiblement a 1'equilibre a chaque instant et que 
la composition du melange polyphasique est variable tout le long de la conduite, la 
masse de chaque constituant du melange etant definie globalement par une equation 
5 de conservation de la masse sans consideration de son etat de phase, et en ce que Ton 
utilise par exemple un schema numerique explicite en temps de fa$on a faciliter la 
resolution des equations du modele. 

Le traitement de l'apparition et la disparition des phases est rendu plus simple 
par cette approche compositionnelle car la masse de chaque constituant y est 
10 consideree globalement sans tenir compte de son etat de phase (monophasique ou 
polyphasique). On evite ainsi les difficultes liees a l'approche "par phase" precedente 
ou il y a une equation de conservation par phase et ou done le nombre d'equations de 
conservation de masse varie a chaque changement d'etat des constituants selon 
qu'une phase apparait ou disparait. 

15 La methode proposee permet de traiter facilement les phenomenes 

d'apparition et de disparition de phases sans rencontrer les problemes de convergence 
des solutions qui se produisent parfois avec les methodes existantes ce qui assure la 
robustesse du code. 

Suivant un mode de mise en ceuvre, la methode comporte une resolution des 
20 equations de transfert d'energie decouplee de celle concernant la conservation de la 
masse et la quantite de mouvement. 

La methode permet une resolution des equations de transfert d'energie, 
decouplee de celle concernant la conservation de la masse et la quantite de 
mouvement et elle comporte de preference I'utilisation d r un schema numerique 

25 explicite en temps, avec cette consequence tres utile que les masses de chacun des 
constituants sont le resultat du schema numerique sans aucun appel iteratif au 
modele hydrodynamique. Connaissant les masses des constituants et la temperature, 
le modele thermodynamique integre determine la pression et la composition du 
melange et notamment la fraction volumique des phases. La detection de Tapparition 

30 et de la disparition des phases est ainsi rendue plus robuste. 
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Avec ce decouplage permis par la methode, on evite ainsi d'avoir a resoudre 
simultanement le modele thermodynamique et le modele hydrodynamique. Par la, on 
evite les conflits qui peuvent naitre du fait que les solutions que fournissent 
respectivement ces modeles sont a priori d'egale pertinence, et difficiles a accorder 
les unes avec les autres. II en resulte que la detection des phenomenes d'apparition et 
de disparition des phases est simple et robuste. 

Suivant un mode de mise en oeuvre, la methode comporte l'assimilation d'un 
melange multi-composants tel qu'un fluide petrolier s'ecoulant dans des conduites, a 
,un_melange,c£mp£rtant,u^ de composants et par exemple a un 



10 melange binaire equivalent (a deux constituants) ayant sensiblement la meme 
enveloppe de phase que le melange reel, de fagon a rendre moins couteux 
1'etablissement du modele compositionnel. 

La methode comporte avantageusement ['utilisation d'un module integre de 
determination des parametres thermodynamiques (equilibre de phase et proprietes de 
15 transport) qui donne des resultats plus representatifs que ceux tires de tables 
precalculees. 

D'autres caracteristiques et avantages de la methode selon l'invention, 
apparaitront a la lecture de la description ci-apres de modes de realisation decrits a 
titre d'exemples non limitatifs, en se referant aux dessins annexes ou : 

20 - la Fig.l presente les conditions aux limites prises en compte; 

- la Fig. 2 illustre le mode de calcul stationnaire approche dans le cas ou la 
temperature est connue; 

- la Fig.3 presente l'algorithme general du calcul stationnaire dans le cas ou le calcul 
de la temperature est demande; 

25 . — la Fig.4 presente l'algorithme de determination de la pression connaissant les 
masses des constituants en utilisant un flash & pression et temperatures imposees; 



- la Fig.5 montre l'algorithme de determination de la pression connaissant les 
masses des constituants en utilisant un flash a volumes et temperatures imposes; 



6 
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- la Fig. 6 presente 1'algorithme de calcul permettant de determiner les grandeurs 
caracterisant I'ecoulement a partir des grandeurs conservatives fournies par le 
schema numerique lors du calcul transitoire; et 

- la Fig.7 presente 1'algorithme de resolution de I'equation de transferts thermiques. 

5 1) Inconnues et equations 

La realisation du modele n composants, p phases, comporte la resolution 
d'equations de conservation de la masse pour chacun des constituants, de 
conservation de la quantiteUe mouvement du melange et deTenergie du melange, 
que l'on va definir ci-apres, en notant comme suit les differents parametres : 



10 




1.1) Inconnues: 






fraction massique du composant i dans la phase j 




c : 
i 


fraction massique totale du composant i 




*r 


fraction volumique de la phase j 






vitesse de la phase j (m/s) 


15 


P : 


pression (Pa) 




T : 


temperature (K) 






enthalpie massique de la phase j 




Pr 


masse volumique de la phase j (kg/ m 3 ) 




T w : 


frottement a laparoi (Pafm) 


20 


Qw- 


terme d'echange de chaleur a la paroi (w/ m ) 




6 : 


angle de la conduite avec Thorizontale 



V M = — — : vitesse barycentrique du melange (m/s) 

p = Zp ; Rj masse volumique moyenne du melange, (kg/ m 3 ) 
PA 

Xj = - fraction massique de la phase j 

25 g : acceleration de la pesanteur (mj V) 

S : surface de passage du fluide (m 2 ^ 

W : variables conservatives 
F : flux du schema numerique 
Q : termes sources 
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Dans Tapproche compositionnelle selon I'invention, la conservation de la 
masse est verifiee pour chacun des constituants. Le transfert de masse entre phases 
n'apparait pas explicitement dans ces equations mais il intervient par le fait que Ton 
decrit le fluide comrne un melange de composition variable le long de la conduite. 
5 Le melange est suppose etre a Tequilibre a chaque instant. 

1.2) Definitions : 

On d6signera dans la suite de la description de fagon abregee par : 

- "flash" un sous-programme integre de calcul des proprietes thermodynamiq ues 
(equilibre liquide-vapeur, composition de chacune des phases) a l'aide d'une 

10 equation d'etat; 

- "flash (P, T)", un "flash" realise en connaissant la composition globale du 
melange, la pression et la temperature; 

- "flash (T, V)", un "flash" realise en connaissant la composition globale du 
melange, la temperature et la masse de chacun des constituants avec une 

15 determination de la pression au cours du calcul, de fagon que les masses soient 



- "flash (P, H)", un "flash" realise en connaissant la composition globale du 
melange, la pression et Tenthalpie totale, avec une determination de la temperature 
au cours du calcul, de fagon que l'enthalpie soit verifiee. 



verifiees; 



20 



1.3) Equations 



Les equations de conservation traitees sont les suivantes : 



- une equation de conservation de la masse pour chaque constituant i : 




- une equation de conservation de quantite de mouvement du melange : 



25 




- une equation d'energie du melange : 



8 
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ou Q,, est le terme de flux d'energie a la paroi. 

Pour tenir compte de l'approche compositionnelle choisie, on impose en plus 

que : 

5 t x ) = 1 V ; 

i' - i 

p . 
c, =— 

I =1 

7=1 

Comme conditions aux limites, on peut imposer par exemple a l'amont les 
debits massiques de chaque constituant et la temperature, et a 1'aval, la pression. 

10 Coraportement thermodynamique du melange 

Par application des lois de la thermodynamique, on obtient les proprietes 
physiques du fluide necessaires au code compositionnel. 

Pour une pression, une temperature et une composition globale donnees, la 
thermodynamique permet de connaitre les fractions massiques de chacun des 

15 composants dans chacune des phases. En supplement, elle calcule la masse 
volumique de chacune des phases presentes et elle permet de deduire les taux 
volumiques de chacune des phases dans le melange global et ainsi de detecter si le 
melange est diphasique, monophasique liquide ou monophasique gaz. Des lois 
elementaires qui seront definies plus loin permettent le calcul des proprietes de 

20 transport de chacune des phases : viscosite, conductivity thermique, chaleur 
specifique, enthalpie massique, ainsi que la tension interfaciale. 



9 
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La loi thermodynamique globale (calcul de I'equilibre et des proprietes du 
melange) est notee : 



i = l,n 
J J = hp 



proprietes de transport 
K proprietes thermiques 

Loi hydrodynaimique 



5 Dans le cas d'un melange, le comportement hydrodynamique est caracterise 

par Timportance du glissement entre les phases en presence : c'est-a-dire par I'ecart 
entre les vitesses de la phase gaz et de la phase liquide dans le cas d'un melange 
diphasique. Le glissement depend des proprietes thermodynamiques des fluides, de 
la fraction massique de gaz, de la vitesse moyenne du melange. II est calcule par une 
10 fonction appelee fonction hydrodynamique qui determine egalement la configuration 
de l'ecoulement et les termes de frottement. Cette fonction est notee : 



ou dV v =V,-V, 

Pour constituer le modele physique convenant a un ecoulement diphasique, 
15 on utilise par exemple les lois de fermetures suivantes bien connues des specialistes : 

- pour le frottement parietal on utilise les coefficients de frottements de type 
Churchill en ecoulement turbulent et de type Poiseuille en ecoulement laminaire ; 

- pour le frottement interfacial, on utilise une loi similaire a celle de Andritsos, N. et 
Hanratty T.J.,1987, Influence of interfacial waves in stratified gas-liquid flows, 

20 AiChe Vol. 33, pp 444 - 454 ; 

- pour le diametre des bulles, on utilise une loi du type de celle proposee par Hinze, 
J.O, 1955, Fundamentals of the hydrodynamic mechanism of splitting in 
dispersion processes, AiChe J. f Voll , pp 289-295 ; 



10 
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- pour la fraction volumique de gaz dans les bouchons liquide, on utilise une loi 
inspiree de Andreussi P. , Bendiksen K., 1989, An investigation of void fraction in 
liquid slugs for horizontal and inclined gas-liquid flow, Int. J. Multiphase Flow, 
Vol. 15-2, pp 937-946. 


5 


Module thermique 

Le terme Q w dar\s l'equation des transferts thermiques correspond au terme 
d'echange du a la contribution des differents modes de transfert thermiques: la 
conduction a travers la conduite (la paroi et les isolants), la convection au sein du 


10 


fluide et 1'echange entre le fluide et le milieu environnant (air, sol ou rner). Pour la 
convection au sein du fluide, on tient compte de la configuration de recoupment. 

II) Proprietes physiques des fluides 

Pour caracteriser les melanges, on classe leurs grandeurs caracteristiques en 
trois groupes : 

o Les informations resultant de 1'etude de 1'equilibre entre les phases: 


15 


o L'etat du melange (p-phasique, triphasique, diphasique, 
monophasique gaz ou monophasique liquide); 

o Les fractions massiques de chacun des composants dans chacune 

des phases et le taux volumique de chaque phase dans le melange; 
o Les masses volumiques de chacune des phases. 


20 


o Les proprietes de transport (utiles pour la resolution hydrodynamique): 

o Les viscosites de chacune des phases; 
o La tension interfaciale. 


25 


o Les proprietes utiles pour la moderation des transferts de chaleur: 

o Les enthalpies massiques de chaque phase; 

o Les conductivit6s thfermiques dans chaque phase; 

o Les chaleurs specifiques des phases. 

Comme outils pour la modelisation thermodynamique, on peut utiliser : 
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a) des correlations: c'est-a-dire des lois simples permettant de representer 
quantitativement les phenomenes physiques dont l'avantage est d'etre tres peu 
couteuses en temps de calcul, ou bien 

b) des proprietes calculees a partir d'un programme thermodynamique 
5 complet avec resolution d'une equation d'etat: soit par 1'intermediaire d'une table 

remplie par un traitement prealable, soit par 1'appel a un "flash" integre a chaque fois 
que Ton a besoin de connaitre les caracteristiques du fluide. 

Presentation des correlations 

Les proprietes physiques : masses volumiques p, viscosites, et les tensions 
10 interfaciales a sont calculees par les formules algebriques simples. Par exemple dans 
le cas d'un melange diphasique gaz-liquide, le gaz aura un comportement proche de 
celui d'un gaz parfait et Ton utilisera les relations suivantes : 

PG ~ PGNorm " p T 
r Norm ' i 

V G = V GNorm 

a L 

V L = V LNorm 

Les valeurs p CNom , P Nom , T Nom ,p LNom , V GNorm , V LNom ,a u <J Nonn sont 

foumies par I'utilisateur, pour chaque simulation, de fa$on a reproduire au mieux le 
20 comportement du fluide modelise. 

Coefficients d'equilibre 



On ne modelise pas directement les fractions massiques de liquide et de gaz 

P°i 

de chacun des constituants mais les facteurs d'equilibre K h K t = — 1 pour chacun 
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des constituants, ou P°i est la pression de saturation a temperature donnee du 
composant i. 

Apparition et disparition des phases : 

Pour determiner 1'etat du melange, deux notions interviennent: la pression de 
5 bulle et la pression de rosee. A composition donn6e, les pressions de bulle et de 
rosee se calculent facilement. L'etat du melange depend alors de la valeur de la 
pression: 

te-m61ange^est-diphasique^si P^ < P-<-P bTltre 

le melange est monophasique gaz si P < P msie 

1 0 le melange est monophasique liquide si P > P huUe 

Calcul des proprietes a Tissue du flash 

Pour le calcul des proprietes physiques, on utilise par exemple le modele 
thermodynamique de Peng-Robinson avec translation de volume, bien connu des 
specialistes. Les viscosites sont calculees par la methode dite de Lohrentz Bray 

15 Clarck, les chaleurs rnassiques a pression constante et enthalpies, par la methode dite 
de Passut et Danner, qui est de forme polynomiale en fonction de la temperature, et 
utilise sept coefficients, la tension interfaciale, par la methode dite des parachors, ces 
parachors et l'exposant caracteristique etant calcules par la methode dite de Broseta, 
toutes methodes qui sont egalement bien connues des specialistes et decrites par 

20 exemple dans les publications suivantes : 

^ Broseta D., et al., 1995, Parachors in Term of Critical Temperature, Critical 

Pressure and Acentric Factor, SPE Annual Technical Conference and Exhibition, 

Dallas, USA, 22-25 oct 1995 ; 
- Passut C.A., et aL, 1972, 1& EC, Process des, dev., J J, 543 {1972) ; 
25 - Peneloux A., et al, 1982, A consistent correction for Redlich-Kwong-Soave 

volumes, Fluid Phase Equilibria t 8 (1982) pp 7J23 Elseviers Science Publishers 

(Amsterdam) ; 
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- Peng D.Y., et al., 1976, A new two-constant equation of state. Ind. Eng. Chem. 
Fund. 15, 59-64 (1976). 

Assimilation a un melange a nombre reduit de composants 

Comme on Ta vu, une caracteristique importante de la methode selon 
5 ['invention est la possibility quelle offre d'assimiler les melanges multi-constituants 
tels que les fluides petroliers par exemple, a un melange d'un nombre plus reduit de 
pseudo-constituants dont les proprietes sont les plus proches possibles de celles du 

metei^e-reely^eei-a^ar^ 

complexe, et par exemple a un melange binaire de deux pseudo-composants. 

] 0 Utilisation des tables de proprietes 

On utilise de preference des "flash" integres pour ce qui concerne notamment 
la determination de la fraction massique vapeur, et surtout si Ton veut traiter des 
> fluides a plus de deux constituants, qui donnent des resultats bien plus representatifs 
que des tables de proprietes thermodynamiques pre-remplies a Taide d'un programme 
15 de calcul sur la base de la description binaire du fluide. 

IE) Determination approchee de Tetat stationnaire 

Avant toute simulation, on procede a un calcul approche de I'etat initial 
stationnaire. Ce calcul permet, en partant d'une solution proche de I'etat initial reel, 
de reduire le temps de convergence pour atteindre un etat stationnaire en accord avec 
20 le schema numerique. 

— " ■ Pour obtenir cet etat, les equations sont resolues en ne considerant pas les 
termes en derivees temporelles et en negligeant les termes d T inertie dans l'equation de 
quantite de mouvement Les donnees utilisees sont les conditions aux limites 
indiquees sur la Fig.l, ou T est la temperature, ql,..., qn les debits massiques des 
25 constituants et Pest la pression. 

On resout le systeme d'equations suivant : 
p 

X PjRjVjX^q' /S pouri=l...n 
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- = r w -(lp,^)^sin0 

I IPjRjHjVj =1(RjVj)o^ -Q„ +V m T„ 

^(v M , Xiy r(p,T,c it ...,c n _ p+l ),dv) = o 





proprietes de transport 






^ proprietes thenniques J 





5 Les inconnues de ce probleme ainsi pose sont : P y T,c ] ,...,c„_ /;+ , ,dV u ,Vj . 



La demarche suivie consiste a effectuer un premier calcul de l'aval vers 
l'amont, qui permet de determiner la pression a Famont pour un profil de temperature 
impose. Seuls les calculs de 1'hydrodynamique et de la thermodynamique sont 
effectues. En partant de l'aval (Fig.2), la pression est imposee, les debits sont connus 
10 (egaux a ceux de I'aniont), la temperature est imposee (cas de profil impose, ou 
estimee dans le cas contraire). 

Si le profil de temperature doit etre calcule, l'algorithme de calcul, est : 

o Calcul de l'aval vers l'amont en ayant estime une temperature a l'aval 
en fonetion des donnees sur Tenvironnement de la conduite. 

15 o Calculs de l'amont vers l'aval effectues en estimant la pression amont 

et en effectuant les calculs thermodynamique, hydrodynarnique et 
thermique. La resolution se fait par une methode de calcul approche 
de type Newton a partir de la pression amont afin de retrouver la 
pression imposee a l'aval, comme schematise sur la Fig. 3. 

20 IV) Determination des etats non-stationmaires 

Les comportements instationnaires sont provoques soit par des variations des 
conditions aux limites par rapport & un etat initial stationnaire, soit encore par la 
geometrie irreguliere du terrain ou de Installation qui provoque des ecoulements 
instables que les specialistes designent par "terrain slugging" et "severe slugging". 
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D6couplage de liquation de transferts thermiques : 

Comme la temperature varie beaucoup moins que la pression, la composition 
et les grandeurs hydrodynamiques, il est done possible de resoudre les echanges 
thermiques de fagon decouplee du calcul de la conservation de la masse et de la 
5 quantite de mouvement. 

Le systeme (quantite de mouvement, masse) peut etre resolu par un schema 
numerique connu ecrit sous forme de volumes finis, explicite en temps par exemple, 
d'ordre 1 ou 2 en espace, tel que decrit par Roe P.L., 1980, "The use of Riemann 
problem in finite difference scheme", in Lectures notes of Physics 141. 

10 Liquation des transferts thermiques est resolue par une methode dite des 

caracteristiques. Le decouplage entre la resolution thermique et hydrodynamique, a 
pour consequence qu'a ce stade de la resolution, la temperature est connue. 



On considere les variables conservatives suivantes : 



15 



Les valeurs de ces grandeurs sont fournies par le schema numerique de Roe 
ci-dessus defini. 



Le systeme d'equations a resoudre est le suivant : 



W Mvl = S pjRjVj 



* propriete de transport = T lherm0 \P y T,c x 
^ proprietes thermiques > 




<t>(v M ,x j ,r lhrrm „,dv ij ) = o 
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Resolution pour les bords internes du maillage : 

Le schema numerique utilise permet de decoupler la resolution 
hydrodynamique de la resolution thermodynamique. Le probleme se ramene a la 
determination des grandeurs physiques a partir des grandeurs conservatives fournies 
5 par le schema numerique par la methode decrite ci-apres. 

La connaissance des masses de chacun des constituants (W;) permet de 
determiner la concentration massique de chaque composant i : 

_ ^ , 

^ =~ 

Calcul de la pression et des proprietes thermodynamiques : 

10 On peut utiliser ici une methode standard consistant a iterer sur la pression en 

faisant des appels a un "flash (P, T) M jusqu'a ce que la pression conduise a des masses 
de chacun des constituants egales a celles que le schema numerique a fournies, 
suivant I'organigramme de la Fig.4. 

On peut noter qu'un gain en temps calcul important a ete obtenu dans le cas 
15 d'un appel a un "flash (T, V)" integre, en imposant des valeurs de temperatures et de 
volumes (cf. organigramme de la Fig.5); les iterations sont faites au sein de la 
thermodynamique et evitent des calculs redondants. 

Appel a la fonction hydrodynamique 

Connaissant la pression et toutes les caracteristiques du fluide, il est possible 

W 

20 d'effectuer le calcul de la vitesse barycentrique V M = — ^.'Un appel de la 

P 

fonction hydrodynamique ®{y Ml X j ,r Thtrmo ,dV ij 'j=0, permet de calculer les 
vitesses des phases. 

II est interessant de souligner que la resolution du modele hydrodynamique 
est identique dans la partie transitoire et dans la partie stationnaire. Le diagramme 



I 
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general de resolution permettant d'obtenir les grandeurs physiques a partir des 
grandeurs conservatives est montre a la Fig. 6. 

Traitement des conditions aux limites 

Dans le cas general, il y a n valeurs propres positives et une valeur propre 
5 negative: 

k x < 0 < A 2 < A 3 <...< X n 

VoMifiofTa'l^limite aval 

Dans le cas general, il y a done une caracteristique entrante et n 
caracteristiques sortantes. La condition a la limite traduit rinformation entrante. La 
1 0 pression est imposee 

N equations de compatibility associees aux valeurs propres positives, 
traduisent reformation sortante : 

Pour j = l,n 

15 n 

. IX^+«a,„W Mw =Aj 

La loi hydrodynamique et la loi thermodynamique doivent etre egalement 
verifiees. 

La resolution du systeme non lin6aire ainsi obtenu permet de determiner les 
masses (Wj) et la quantite de mouvement (W^vO- 

20 Condition a la limite amont : 

Dans le cas general, il y a n caracteristiques entrantes et 1 caracteristique 
sortante. Les conditions a la limite, sur les debits imposes, q^q^ ^ qn de chacun 

des composants, traduisent rinformation entrante : 

Les equations a verifier sont done les suivantes : 
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q,.-(T.PjXjVjx'j)-S = 0 pourj=1an 

L'information sortante est traduite par I'equation de compatibilite associee a 
la valeur propre negative : 

Zap, ( 3 ) 

5 La loi hydrodynamique et la loi thermodynamique doivent etre egalement 

verifiees. 

La resolution du systeme non lineaire ainsi obtenu permet de determiner les 
masses (Wj) et la quantite de mouvement (W^vt)- 

Remarque generale : 

10 Les differentes etapes de resolution numerique tant sur les bords que sur les 

conditions aux limites font appel a des calculs de derivees partielles. Par souci de 
robustesse, de precision et de gain en temps calcul, la plupart des derivees sont 
calculees de fagon analytique notamment les derivees des grandeurs 
thermodynamiques. 

1 5 Transferts tthermiques transitoires 

A ce stade de la resolution, la composition du melange et la pression sont 
connues (resolution de la conservation de la masse et de la quantite de mouvement). 

Une methode de caracteristiques (Fig.7) resout Tequation des transferts 
thermiques ; elle fournit la valeur de l'enthalpie massique du melange. 

20 Uappel a la loi thermodynamique permet de definir la temperature : 

o La methode standard consiste a faire des appels successifs a un "flash (P, T) M 
en faisant evoluer la temperature jusqu'a ce que l'enthalpie du melange 
calculee soit identiqxie a celle fournie par le schema numerique. 



25 



o Un gain de temps de calcul peut etre obtenu en ecrivant directement un 
"flash (P, H)". Les iterations, dans ce cas, sont faites au sein meme du 
modele thermodynamique. Le modele thermique fournit le terme d'echange 
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therrniqueQnt • Q u i est ' une des parties du terme source, Q e nth I'equation 
de la chaleur. 

La methode selon I'invention qui permet done de modeliser la variation de la 
composition d'un melange a nombre reduit de composants dans l'espace et au cours 
du temps, a ete validee experimentalement sur des cas reels. 
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REVINDICATIONS 



1) Methode pour constituer un modele representatif des ecoulements 
permanents et transitoires dans des conduites d'un melange polyphasique, comportant 

5 1'utilisation d'un modele hydrodynamique de type a flux de derive et d'un modele 
thermodynamique integre pour definir les proprietes des constituants et la resolution 
d'un ensemble d'equations de conservation de la masse, de conservation de la 

quaiitfte-deTnouvem 

que Ton forme le modele en considerant que le melange est sensiblement a Fequilibre 

10 a chaque instant et que la composition du melange polyphasique est variable tout le 
long de la conduite, la masse de chaque constituant du melange etant definie 
globalement par une equation de conservation de la masse sans consideration de son 
etat de phase. 

2) Methode selon la revendication 1, caracterisee en ce qu'elle comporte 
15 1'utilisation d'un schema numerique explicite en temps de fa$on a faciliter la 

resolution des equations du modele. 

3) Methode selon Tune des revendications 1 ou 2, caracterise en ce qu'elle 
comporte l'assimilation de melanges multi-composants a des melanges constitues 
d'un nombre plus reduit de composants equivalents. 

20 4) Methode selon la revendication precedente, caracterisee en ce qu'elle 

comporte l'assimilation de melanges multi-composants a des melanges binaires 
equivalents. 

5) M6thode selon l'une des revendications precedentes, caracterisee en ce 
qu'elle comporte une resolution des equations de transfert d'energie decouplee et 

25 celle concernant la conservation de la masse et la quantite de rnouvement. 

6) Methode selon l'une des revendications precedentes, caracterisee en ce 
qu'elle comporte 1'utilisation d'un module integre et optimise de determination des 
parametres thermodynamiques definissant l'equilibre des phases et les proprietes de 
transport du melange. 
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